Escape from a zero current state in a one dimensional array of Josephson junctions 
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A long one dimensional array of small Josephson junctions exhibits Coulomb blockade of Cooper 
pair tunneling. This zero current state exists up to a switching voltage, V sw , where there is a sudden 
onset of current. In this paper we present histograms showing how V aw changes with temperature for 
a long array and calculations of the corresponding escape rates. Our analysis of the problem is based 
on the existence of a voltage dependent energy barrier and we do not make any assumptions about 
its shape. The data divides up into two temperature regimes, the higher of which can be explained 
with Kramers thermal escape model. At low temperatures the escape becomes independent of 
temperature. 
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Escape from a metastable state can be used to describe 
a wide variety of physical phenomenon, from solid state 
diffusion jj] to magnetic switching Q , and has for a long 
time attracted attention from both a theoretical and an 
experimental point of view (for a review see e.g. ||). 
Thermal escape over an energy barrier was studied the- 
oretically by Kramers in the 1940s (J]. More recently, 
the influence of quantum tunneling on the escape pro- 
cess has been an important means of studying the quan- 
tum behavior of macroscopic variables. In this context, 
escape from the zero voltage state of a Josephson junc- 
tion has been thoroughly studied |^|, || . In this paper we 
adopt the Kramers model to analyze the escape from the 
zero current state of a one dimensional array of Joseph- 
son junctions exhibiting a Coulomb blockade of Cooper 
pair tunneling. We find that the escape process is well 
described by thermal activation at higher temperatures, 
but at lower temperatures the escape process becomes 
independent of temperature. 

In the Kramers thermal activation model, the rate of 
escape F, over an energy barrier At/, is given by, 



oj p , At/ 

k — — exp 

2tt PV k B T' 



(1) 



where ks is Boltzmann's constant and uj p is the attempt 
frequency, or the frequency of small oscillations around 
the metastable minimum. The dimensionless transmis- 
sion factor k depends on the shape of the barrier and 
the damping M. This model has been applied to single 
Josephson junctions, where the energy barrier is a one 
dimensional tilted washboard potential, and the tilt is 
controlled by bias current applied to the junction, AU(I) 
[& 0, Mj. Long Josephson junctions (modeled by parallel 
arrays of single junctions) have also been studied, and in 
this case the escape is from a metastable minimum in a 
multi-dimensional space 0. All Josephson junction sys- 
tems studied with this technique thus far exhibit a zero 
voltage state with a finite supercurrent, above which the 
system switches to a finite voltage state. A voltage across 



the junction develops when the superconducting phase, 
4>, across the junction escapes over the Josephson po- 
tential energy barrier to a free-running state (V oc cf>). 
The escape can be due to thermal activation or quan- 
tum tunneling of the junction phase through the energy 
barrier. The important quantities for the escape problem 
are the height of the Josephson potential and the attempt 
frequency, which are given in terms of other junction pa- 
rameters such as normal state resistance, superconduct- 
ing energy gap and junction capacitance. This escape 
description applied to the Josephson junction is based 
on the classical dynamics of the Josephson phase, as em- 
bodied in the Josephson effect. 

In this letter we analyze the Coulomb blockade of 
Cooper pair tunneling, a phenomena which is in many 
ways dual to the Josephson effect || [l(J. In one- 
dimensional arrays of small capacitance Josephson junc- 
tions, the coupling of the Josephson phase to dissipation 
can be manipulated so that there are very large quantum 
fluctuations of the phase 



11 . In this case, the number 



of Cooper pairs becomes a classical variable, and trans- 
port through the array is described by a Coulomb block- 
ade, where a zero current state persists until a switching 
voltage is exceeded. One- and two-dimensional arrays 
of small capacitance Josephson junctions have been well 
studied in the context of the superconductor-insulator 
transition Here we focus on analysis of transport 

properties of a ID array deep in the insulating regime. 
We do this by studying fluctuations of the switching volt- 
age and analyizing them in the context of the Kramers 
model. This is the first time to our knowledge that such 
an experimental study has been carried out on Cooper 
pair tunneling in the insulating regime. 

The array is made of Al, with AI2O3 tunnel barriers, 
on a Si02 surface using the standard two angle evapora- 
tion technique [ p"3| . Electron beam lithography was used 
to define the mask pattern, which made 400 junctions 
in series. Each junction had dimensions lOOnmx 150nm, 
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and each island was connected by two junctions in par- 
allel, forming a Superconducting QUantum Interference 
Device (SQUID). The charging energy E c of the SQUID 
was Ec = e 2 /2C = 59/ieV, where we have assumed a spe- 
cific capacitance of 45 jFj fim 2 |l4| . The junctions normal 
state resistance R n controls the Josephson energy at zero 
magnetic field, E J0 = [Rq/Ru) A/2 where A = 200/zeV 
is the superconducting energy gap of Al. For the ar- 
ray presented here we had R n = 4.6k£l per SQUID, or 
Ejo = 142/ieV. 

The advantage of the SQUID geometry is that we can 
tune the Josephson coupling by applying a magnetic flux 
to the SQUID loops Ej = E JQ \ cos(ttBA/§ ) \, where B 
is the magnetic field perpendicular to the SQUID loop, A 
is the effective area of the loop (in this case 0.85 x 0.2/im 2 ) 
and <&o = h/2e = 2.06 x 10~ 15 Wb is the superconducting 
flux quantum. 

The arrays were measured in a dilution refrigera- 
tor with well filtered measurement lines. The chip is 
mounted on a socket which is in turn mounted in a RF 
tight copper can. The leads into the can are 1 meter 
of Thermocoax which attenuates at high frequency [ p"5| . 
The can and filtered leads are mounted at the lowest 
temperature. The bonding pads and contact leads on 
the chip were fabricated of Au to within 50/im of the ar- 
ray edge. Low noise and high input impedance amplifiers 
were used to measure the voltage and current. 

In Fig. |l| we see a typical current-voltage (I-V) char- 
acteristic of a ID array in the Coulomb blockade state. 
When increasing the bias voltage across the array the 
Coulomb blockade state is interrupted by a sudden onset 
of current at the switching voltage, V sw . Subsequently 
lowering the bias voltage across the array, the Cooper 
pairs are retrapped at a voltage V r < V sw , giving rise to 
a hysteresis in the current-voltage characteristics. The 
retrapping voltage remains the same at fixed magnetic 
field, whereas the switching voltage fluctuates. Both V sw 
and V r depend on Ej and are modulated with magnetic 
field having a period associated with one flux quantum 
in each SQUID loop, as can be seen in the inset of Fig. 
[j]. This modulation is clear evidence that the switching 
voltage is associated with the Cooper pair tunneling in 
the array. 

By making 10 4 sweeps of the bias voltage and each 
time registering the switching voltage V sw , we obtain a 
switching histogram representing the probability distri- 
bution P(V) from which the escape rate can be calculated 
through the formula H, 



T(V) 



1 dV 

~Kv~dT 



in 



\v>V v>V+AV 



P(v) (2) 



where dV/dt is the ramp rate of the voltage and AU is 
the width of the bins in the histogram. 

To establish that the switching process is indeed sen- 
sitive to the applied voltage, we have varied the sweep 



rate by changing the frequency of our triangular sweep 
between 5 and 200Hz, while keeping the amplitude of 
the sweep constant at 1.2mV. For lower sweep rates the 
histogram shifted to lower voltages as we are less likely 
to reach higher voltages before escape occurs. However, 
when the escape rate T(V) was the calculated according 
to (Q) and comparison was made for different sweep rates, 
we found see that for sweep rates lower than 0.12V/s, the 
escape rate at a given voltage was essentially unaffected 
by the sweep rate. We also found that for sweep rates 
higher than 0.12V/s, the escape rate T(V) becomes sen- 
sitive to the sweep rate. This limitation on sweep rate 
was set by the biasing and current measurement circuit 
together with the stray capacitance of the cryostat leads. 
From this analysis we conclude that the switching pro- 
cess, when measured with sweep speed below 0.12V/s, is 
truly a voltage dependent escape problem. We also in- 
vestigated if the switching process is affected by the time 
spent in the dissipative (finite current) state by putting 
in a delay and waiting at zero voltage between succes- 
sive ramps. The result of this investigation was that the 
switching histogram for fixed sweep rate and temperature 
remained the same, independent of delay time. 

Figure || shows a selected number of histograms at dif- 
ferent temperatures. The inset of Fig. || shows how the 
voltage value at the peak of the histogram changes ver- 
sus temperature. At temperatures above about 130mK, 
the displacement of the peak position is proportional to 
the temperature. We also see that the width of the dis- 
tribution increases as the temperature is decreased. The 
retrapping voltage put a lower limit on the switching volt- 
age and thereby sets a maximum temperature of about 
180mK. Therefore we limit ourselves to temperatures be- 
tween 170mK and the base temperature of the cryostat, 
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FIG. 1: Current voltage characteristics of a typical long one 
dimensional array of Josephson junctions. Voltage bias gives a 
hysteretic I — V curve, with a switching voltage, Vsw, and the 
retrapping voltage, V r . The arrows on the curve indicate the 
direction of the sweep. The inset shows how V r is modulated 
with B. 
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FIG. 2: A selection of histograms containing 10 4 switching 
events, each at different temperatures. As the temperature 
increases the histogram shifts towards lower voltages. At high 
temperatures it is evident how the retrapping current affects 
the histogram as V sw gets closer to V T . The inset shows how 
the peak of the histogram moves with temperature. 
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FIG. 3: lnF versus the inverse temperature at I.fmV. Two 
temperature regimes are seen, the higher of which follows the 
thermal escape model. 



about 50mK. 

From the histograms we can determine the escape 
rate r(V) via equation (||). In Fig. || we have plotted 
T(V = l.lmV) versus 1/T. Here we see that in the 
high temperature region between 170mK and 130mK, the 
data is well described by Kramers' model of thermally 
activited escape, equation (|l]). However at low tempera- 
tures, below ca 130mK, the escape rate crosses over to a 
different behavior, becoming independent of temperature 
at the lowest temperatue attainable. By extrapolating 
the linear fit in the low temperature region shown in Fig. 
|3|, we find the intersection with the In T-axis and make an 
estimate of the factor, kuj p of the order of 10 9 — 10 10 Hz. 
An alternative way to find nu> pi which uses the data at 
all voltages, is to extract the energy barrier AU(V) from 
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FIG. 4: Selected curves showing how all data can be reduced 
to a single AU(V) curve in eq. ([j]), with ku p — 80GHz. The 
low temperature curves with the open symbols are fitted with 
an effective temperatures indicated in the parenthesis. 



the escape rate r(V) via equation (|l]). We then adjust 
Kujp so that the data at all temperatures fall on to a sin- 
gle AU(V) curve. We find that in the high temperature 
region between 170mK and 130mK, the data collapse on 
to one curve for kuj p = 80GHz as is shown in Fig. [|. 

We note that this method of determining ku> p differs 
from the traditional analysis of escape from the super- 
conducting state of Josephson junctions. In the super- 
conducting case we have a well established model for the 
energy barrier, and to an excellent approximation, escape 
from the cubic potential can be assumed. The cubic po- 
tential results in a straight line plot of (ln(wp/27rr)) 2 / 3 
versus the bias current || . If we follow this analysis and 
plot (ln(w p /27rr)) 2 / 3 versus bias voltage in the tempera- 
ture range 130 — 170mK using the value klo p — 80GHz, 
we also find a straight line. However, the limited width of 
the distribution and the small temperature range make it 
difficult to say if the data follows the theoretical predic- 
tions (slope oc T -2 / 3 and intersect with the x-axis at Vq). 
We can, however, conclude that the qualitative behavior 
is correct. 

Thus the data are well described by the thermal ac- 
tivation model in the high temperature regime with 
Kui p = 80GHz. Fixing the value of klo p — 80GHz, we can 
then analyze the low temperature regime (T < 130mK) 
by adjusting the temperature in equation ([l]) to find an 
effective temperature. Following this procedure, all data 
for all voltage and tempeature can be reduced to one 
AU (V) as shown for a few selected temperatures in Fig. 
[|. The effective temperature, T e //, is plotted versus the 
thermometer temperature in Fig. pi Here we can see how 
the data can be described by the thermal escape model if 
an effective temperature is assumed that reaches a min- 
imum value of 122mK as the cryostat is cooled towards 
base temperature. 
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There can be several reasons for deviation from the 
thermal activation model at low temperatures. The es- 
cape temperature of 122mK may be a result some source 
of heating of the sample, so that the sample is not in 
equilibrium with our thermometer. We can rule out self- 
heating by the current in the sample |l6| because we are 
measuring escape from a zero current state, where there 
is no dissipation in the sample. The low temperature es- 
cape may also be influenced by RF and microwave noise 
reaching the sample through the filtered measurement 
leads. Another possibility is that the escape process at 
low temperatures is determined by macroscopic quantum 
tunneling (MQT). 

MQT rates have been calculated for the cubic potential 
at T = Ol, 



7.2(1 + 0.87/Q) 



(3) 



where Q is the damping factor. Assuming moderate to 
weak damping (Q 3> 1) and k ~ 1 we can calculate 
an MQT escape temperature T esc ~ 85mK, which is in 
reasonable agreement with the effective escape tempera- 
ture found with our method of analysis. However, simi- 
lar measurements and analysis on a slightly shorter array 
found a much lower kuj p « 10MHz. This would result in 
T eS c ~ 10/iK in the MQT model, in poor agreement with 
the effective escape temperature T e ff(0) w lOOmK found 
in the data analysis. 

In conclusion, we have studied the onset of finite cur- 
rent in long, one-dimensional Josephson junction arrays 
in the Coulomb blockade state. By analyzing the temper- 
ature dependence of fluctuations of the switching voltage 
for the onset of finite current, we have shown that this 
onset can be described by the Kramers model for ther- 
mally activated escape over an energy barrier, where the 
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FIG. 5: The effective temperature, T e //, is plotted versus the 
thermometer temperature, T. At high temperatures T e ff = 
T. At low temperatures T e ff was inferred by mapping all 
data to the same AU(V) curve. 



barrier height depends on the bias voltage. The data 
fall into two temperature regimes. At high temperature 
thermal activation applies, but at low temperature an 
effective temperature must be assumed to fit the data 
to the Kramers model. The low temperature behavior 
could be the result of external noise causing the escape 
over the barrier, or possibly, macroscopic quantum tun- 
neling through the barrier. 



We would like to acknowledge helpful discussions with 
R. Kautz and T. Duty, and financial support from the 
NFR, Goran Gustafsson Foundation and the EU project 
SQUBIT. 



[1] 

[2] 

[3. 

[4 
[5 

[6 

[-: 

[8 
[9 

[io; 
in 

[12 
[13 

[14; 

[15, 
[16 

[17] 



Electronic address: karina@nanophys.kth.se 

R. DiFoggio and R. Corner, Phys. Rev. B 25, 3490 

(1982). 

J. Z. Sun, J. C. Slonczewski, P. L. Trouilloud, D. Abra- 
ham, I. Bacchus, W. J. Gallagher, J. Hummel, Y. Lu, 

G. Wright, S. S. P. Parkin, et al., Appl. Phys. Lett. 78, 
4004 (2001). 

P. Hanggi, P.Talkner, and M.Borkovec, Rev. Mod. Phys. 
62, 251 (1990). 

H. A. Kramers, Physica (Utrecht) 7, 284 (1940). 

J. M. Martinis, M. H. Devoret, and J. Clarke, Phys. Rev. 
B 35, 4682 (1987). 

M. H. Devoret, D. Esteve, C. Urbina, J. Martinis, A. Cle- 
land, and J. Clark, Quantum Tunnelling in Condensed 
Matter (Elsevier Science Publishers, 1992), chap. 6, pp. 
313-345. 

T. A. Fulton and L. N. Dunkleberger, Phys. Rev. B 9, 
4760 (1974). 

M. G. Castellano, G. Torrioli, C. Cosmelli, A. Costantini, 
F. Chiarello, P. Carelli, G. Rotoli, M. Cirillo, and R. L. 
Kautz, Phys. Rev. B 54, 15417 (1996). 
P. Agren, K. Andersson, and D. B. Haviland, J. Low 
Temp. Phys. 124, 291 (2001). 

P. Agren, Cooper pair transport and Coulomb blockade 
in one dimensional Josephson junction arrays, Licen- 
tiate thesis. KTH, Stockho lm. Sweden. Available at: 



http://www.nanophys.kth.se (2000) 



D. Haviland, K. Andersson, and P. Agren, J. Low Temp. 
Phys. 118, 733 (2000). 

R. Fazio and H. Van Der Zant, Physics Reports 355 
(2001). 

K. Andersson, A Cooper pair turnstile made from a 
one dimensional array of Josephson junctions, Licen- 
tiate thesis, KTH, Stockholm, Sweden. Available at: 



http://www.nanophys.kth.se (1999) 



P. Delsing, C. D. Chen, D. B. Haviland, Y. Harada, and 

T. Claeson, Phys. Rev. B 50, 3959 (1994). 

A. B. Zorin, Rev. Sci. Instrum. 66, 4296 (1995). 

R. L. Kautz, G. Zimmerli, and J. M. Martinis, J. Appl. 

Phys. 73, 2386 (1993). 

M. H. Devoret, J. M. Martinis, and J. Clarke, Phys. Rev. 
Lett. 55, 1908 (1985). 



